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1. Introduction 



Let (M, g) be a compact Riemannian manifold with boundary. The geodesic ray transform 
/ of symmetric 2-tensor fields / is given by 



where 7 runs over the set of all geodesies with endpoints on dM. All potential fields dv given 
by (dv)ij = !(VjUj + ^ 'j v i) w hh v = on dM belong to the kernel of I. The ray transform 
/ is called s-injective if this is the only obstruction to injectivity, i.e., if If = implies that 
/ is potential. S-injectivity can only hold under certain assumptions on (M,g). A natural 
conjecture is that it holds on simple manifolds, see the definition below. So far it is known 
to be true for some classes of simple manifolds only, including generic simple manifolds, see 



In the cases where s-injectivity is known, there is also a stability estimate that is not 
sharp. In [TT], it is of conditional type with a loss of a derivative, see below. In [16j, 
the estimate is not of conditional type but there is still a loss of a derivative, see ([3]) below. 
On the other hand, if / is a function, or an 1-tensor (an 1-form), there is a sharp estimate, 
see [15] . The purpose of this paper is to prove a sharp estimate for the ray transform of 
2-tensors. 

The geodesic ray transform is a linearization of the boundary distance function and plays 
an important role in the inverse kinematic problem (known also as boundary or lens rigidity), 
see e.g., [TTl [TBI, [Tg| [TT] and the references there. There, one wants to recover (M,g) given 
the distance function on dM x dM or the scattering relation a : (x, £) 1— > (y, rf) that maps a 
given x G dM and a given incident direction £ to the exit point y and the exit direction 77 of 
the geodesic issued from (x,£). 



Definition 1. We say that a compact Riemannian manifold (M,g) with boundary is simple 



(a) The boundary dM is strictly convex, i.e., (V^u, £) > for each £ G T x (dM) where v 
is the unit outward normal to the boundary. 

(b) The map exp x . : exp" 1 M —>■ M is a diffeomorphism for each x G M. 



(1) 
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Condition (b) implies that each pair of points is connected by a unique geodesic depending 
smoothly on the endpoints. It also implies that M is diffeomorphic to a ball, so we can work 
in one fixed chart only. We will fix M fixed, and choose different metrics on it. If (M,g) is 
simple, then we call g simple. 

It is known [IT], see also [IB] , that each symmetric 2-tensor field / G L 2 (M) admits an 
orthogonal decomposition 

f = r + dv, 

where v is 1-form in Hq(M) (vanishing on dM), and / is divergence free, i.e., 5f = 0, where 
(Sf)j = V l /ij, an d V is the covariant derivative. The 1-form v solves 5dv = 8f, v\qm = 0. 
The latter is an elliptic system and the Dirichlet boundary condition is a regular one for it, 
see [UJ [16]. S-injectivity then is equivalent to the following: lf = implies f s = 0. 
Set 

d±SM := {0, u) G TM; x G dM, \u\ = 1, ±(u, u) > 0} , 
where u(x) is the outer unit normal to dM (normal w.r.t. g, of course). Here and in what 
follows, we denote by (u,u) the inner product of the vectors uo, v, and |u;| is meant w.r.t. 
g. Let 7^,0) (£) be the (unit speed) geodesic through (x,uj), defined on its maximal interval 
contained in [0, oo). One can then parametrize all maximal (directed) geodesies in M by 
points on dSM, and with some abuse of notation we denote If{x,u) = If(j x ^>)- 

One of the methods to study s-injectivity of I is the energy estimates method that goes 
back to Mukhometov [El HUE], see also [2], where, for simple manifolds, injectivity of / acting 
on functions / is proved. S-injectivity (injectivity up to d<ft, where = on dM) for 1-forms 
/ is established in [JJ by a modification of the same method. The case of 2-tensors is harder. 
Using the so-called Pestov identity, Sharafutdinov [UJ showed that I is s-injective under an 
explicit a priori bound on the positive part of the curvature of g, that in particular implies 
simplicity of g but it is not equivalent to it. This generalized earlier results on negatively 
curved manifolds [9]. The Pestov-Sharafutdinov approach implies the following stability 
estimate [UJ under the small curvature assumption: 

(2) ||/l! 2(n) < C( \\jJ\ dM \y {dM) \\If\W^sM) + \\If\\ 2 m [a _sM)), V/ G H\M), 

where [j v f)i = fijV 3 . The Sobolev spaces on the r.h.s. are taken with respect to the induced 
measure da(x,ou). In semigeodesic local coordinates x = (x',x n ), the latter is given by 
da(x,u) = (det gfl 2 dx l . . . dx™" 1 (det g) l l 2 da x {uj), where da x (tu) = i| E^-l^V^ 1 A 

• • • A d£> A • • • A d£ n | . 

This estimate is of conditional type — while it implies s-injectivity (under the curvature 
condition), it says that f s is small if If is small and we have an a priori bound on /. 
Moreover, there is a loss of one derivative: the r.h.s. is finite under the condition that 
/ G H 1 , while in the l.h.s., we have only the L 2 norm of f s . 

In [T5 j [T6l 118], the author and G. Uhlmann studied this problem from microlocal point of 
view. Introduce the following measure dfj,(x,u) = \(u,tj)\da(x,uj) on d^SM. It is easy to 
show that / : L 2 (M) — > L 2 (d_SM, dfi) is bounded [JTJ (the first space is a space of tensors, 
actually). The normal operator N = 1*1 is well defined on L 2 (M) then. We extend slightly 
M to a larger manifold with boundary Mi that is still simple so that Mi =>> M. We extend 
tensors defined in M as zero to M\. Then / : L 2 (Mi) — > L 2 (d-SMi, d/i), and one can define 
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1*1 related to Mi that a priori is different from N. On the other hand, it is easy to see that 
when restricted to tensors supported in M, it coincides with N. Then N : L 2 (M) —>■ L 2 (Mi) 
is s-injective if and only if / is s-injective. 

The main results in [151 05] concerning the tensor tomography problem are the following. 
The operator I is s-injective for real analytic simple metrics. Moreover, the set of simple 
C fc metrics, where k 3> 1 is fixed, for which / is s-injective, is open and dense in C k . Therefore, 
we get s-injectivity for a generic set of simple metrics. Next, for any simple g for which I is 
s-injective, one has the stability estimate 

(3) llribcM) < C\\Nf\\ S2{Mi) , V/ G H\M). 

Here H 2 (M) is defined as follows. To the usual if 1 (Mi) norm, we add a term of the kind 

II^V/|| L 2 ([/) + \\x n d x nVf\\ L 2 (u) , 

where (x r , x n ) are semigeodesic local coordinates near dM, and U is any fixed neighborhood 
of dM. The constant C in ([3]) can be chosen locally uniform for g G C k . 

A natural conjecture is that / is injective for all simple metrics. This is still an open 
problem. For any simple metric however, we have an estimate of the kind ([3]) plus the term 
II Kf || l 2 (m) , where if is a smoothing operator. In Proposition [1] below we prove a sharper 
estimate of this kind. It is also known that the solenoidal tensors on the kernel of I form 
a finitely dimensional space of smooth tensors [121 H51 (THl El [TOj . This also follows directly 
from the analysis below, since the inversion problem is reduced to a Fredholm one. 

Estimate ([3]) is not of conditional type anymore but there is still a loss of one derivative. 
Indeed, N is a \l/DO of order —1, and the natural norm on the r.h.s. of ([3]) would be the 
if 1 (Mi) one. Our main result shows that this is the case, indeed. 

Theorem 1. Let g G C k (M), fc > 1, k a simple metric on M, and assume that I is 
s-injective. 

(a) Then 

(4) ||/ - \\lhm)/C < \\Nf\\ m{Ml) < C\\f s \\ LHM) 
with some C > 0. 

(b) The constant C can be chosen uniformly under a small C 3 (M) perturbation of g. 

As pointed out above, the s-injectivity assumption is generically true for simple metrics, 
and holds in particular for metrics close enough to analytic ones [16] or for metrics with an 
explicit bound on the curvature [TT] . 

The new ingredient of the proof is the use of Korn's inequality [191 Collorary 5.12.3], see 

(USD- 

In [T7], the author and G. Uhlmann considered manifolds that are not simple, with possible 
conjugate points, and studied the question of the s-injectivity on / known on a subset T of 
geodesies. The basic assumption is that none of the geodesies in T has conjugate points, 
and the conormal bundles of all 7 G T cover T*M. Under that assumption, results about 
s-injectivity for generic simple metrics, including analytic ones are obtained. A stability 
estimate of the kind ([3]) is also proven there, where iV is modified via a smooth cut-off that 
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restricts the geodesies to T. Without going into detail, we will only mention that Theorem [T] 
generalizes to that case, i.e., the stability estimate in [T7] can be written in the form as 
well. 

Theorem [T] allows us to reduce the smoothness requirement in the generic result in [16J. 

Corollary 1. There exists a dense open set of simple metrics in C 3 (M) so that the corre- 
sponding ray transform I is s-injective (and holds). 

Note that we are not claiming that the set of all C 3 (M) simple metrics with an s-injective 
/ is open. Our success with proving estimate (TjJ that implies the openness depends on 
our ability to show that the problem can be reduced to a Fredholm one. We do this by 
constructing a parametrix, and this requires certain number k of derivatives of g, at least 
k = 2n + 1. Once we have (jl]), we use the singular operator theory to perturb (prj near any 
go G C k with an s-injective J, by C 3 perturbations. 



3. Proofs. 

3.1. Proof of Theorem [T](a). We start with recalling some facts from [TSJ [TB]. We show 
first that 



(NfY> (x) = 2 / w » u> / Ml.UtmUtmjt) dtda x (u), 
Js x m Jo 

where / is supported in M, and we work in M\. Performing a change of variables, we get 
the integral representation 



(5) (Nf) kl (x) 



P(y) dp dp dp dp ^ d 2 (p 2 /2) 

dxdy 



y/detg(x) J p(x, y) n ~ l dy i dyi dx k dx l 



dy, x e Mi, 



where p is the distance function. This form of iV show that iV is a ^DO of order —1 on the 
interior of Mi. Its principal symbol is, see [THl E], 

(6) o- p (N) ijk \x, = 2tt / u^u^d^ ■ u) da x (u), 

Js x M 1 

where ^ ■ uj = C,iW l . This formula generalizes in an obvious way to tensors of any order. It 
follows now easily that iV is elliptic on tensors satisfying Cfij = (solenoidal tensors in the 
Fourier representation), and vanishes on tensors of the type \{^iVj +Cj v i) (potential tensors 
in the Fourier representation). This fact allows us to construct a first order \^DO Q so that 
for any / G L 2 (M), 

(7) QNf = fi h + Kf 

in Mi, where is the solenoidal projection of / (extended as zero outside M) in M 1; and 
K is a compact operator. We can assume that the kernel of Q has a support close enough 
to the diagonal. The need to work in Mi is due to the fact that we can use the (standard) 
^DO calculus in an open set only. For more details, we refer to [T5l [T6] . Note that this 
construction needs only a finitely smooth metric g G C k (M), k ^> 1, that we extend to Mi. 
If we want K to be infinitely smoothing, then we need g G C°°(M). 
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The next step is to construct f s , given fL . This can be done in an explicit way as follows. 
Note that 

(8) f s Mi = Ef s -dw in Mi, 

where E is the extension as zero to Mi \ M, and w solves the elliptic system 

(9) Sdw = SEf s , w\ 9Ml = 0. 

The distribution 5Ef s is supported on dM, and the solution w exists in Hq(Mi), see [16] 
and Lemma [1] below. In particular, w\qm G H l l 2 (dM) is well-defined. If we know w\qm £ 
H 1 / 2 (dM), we can recover w in M because 5dw = in the interior of M, by ([9]). If we 
recover w in M, we recover f s as well, in terms of f^, by (jBJ). Our goal therefore is to 
recover w\sm first. 

We first determine w in Mi \ M, up to a smoothing term, by the relation 

(10) f s Ml = -dw inMi\M, 
see OH]). Since w = on <9Mi, we can integrate the identity 

(11) ^t Wi ^ = [^WlvW 

along geodesies in Mi \ M connecting points on dMi and dM to recover w on dM. Let 
r + (a:,£) > be characterized by 7 X ,^(^) £ <9Mi for t = r + (x,£). Then we get 

= r +a ' C) [/lf 1 ]«(7^(*))7i fC (*)'4€(*)d*. 

for any (x, ^) so that {7 x g()f:), < t < t + (x, ^)} does not intersect M. That also implies 
easily the following non-sharp estimate 

(12) ||HU 2 (M!\A:0 < C||^||i 2 (M) < C||/MjU 2 (Mi\iW)) 

see also [11] for the first inequality. We refer to [T5J [16] for more detail. This approach 
provides also a constructive way to reduce the problem to a Fredholm one. For the proof of 
the theorem however, this is not needed. The new ingredient in this work is that we apply 
Korn's inequality [T^l Collorary 5.12.3], 

(13) IIHIffi(-MiW) < c (||dHU 2 (Mi\M) + IIHU 2 (Mi\m)) • 

This inequality is a consequence of the fact that the Neumann boundary conditions for Sd 
are regular ones. Apply the trace theorem, f|T3|) . ffT2l . and ffTUj) to get 

\\ w \\m/2(dM) < C\\w\\ H i( Ml \M) < C'\\f s Mi \\ L 2( Ml \ M) . 
Now, since w solves the elliptic PDE 5dw = in Mi, we get 

(14) I \W II H 1 (M) < C\ I f Ml 1 1 L 2 (Ml \ A/) , 

see Lemma [1] below. This, together with ([8]) yields, 

(is) ||riU 2 (M) < ll/lf 1 IU 2 (M) + c||/| fl |U WM) < c (||JV/||fli (Ari) + II^/IU^mo) • 

It is worth noting that without the a priori s-injectivity assumption, we got the following. 
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Proposition 1. For any I > 0, there exists k > so that for any simple metric g G C k (M), 

WfWmM) < C {\\Nf\\ HHMl) + \\f\\ H - H M)) , V/ G L 2 (M). 

Estimate (jl5p (or Proposition [1]), together with [T9~| Proposition 5.3.1], implies that if /, 
and therefore, N is s-injective, then there is an estimate as above, with a different C, with 
the last term missing. This completes the proof of Theorem Ufa). 

We return to the elliptic regularity estimate (JHj). If 6d is replaced by the Laplace operator, 
then ( fl4"l) follows from [ 19l Theorem 5.1.3]. If we raise the Sobolev regularity everywhere in 
(117p below by 1, this just follows from the fact that the Dirichlet conditions are regular for 
5b. In our case, we follow the proof of [T9l Theorem 5.1.3] to get the following. 

Lemma 1. Let u G H~ 1 (M), a G if 1 ^ 2 ((5M) be 1-forms. Then the boundary value problem 

(16) Sdw = u in M, w\bm — a 

has a unique solution w G H 1 (M), and the following estimate holds 

(!7) \\w\\ H i (M) < C (||w||h-i ( m) + \W\\m/*(dM)) 

Proof. By a standard argument, first extend a G H l l 2 (M) to a G // 1 (M) by means of a 
fixed bounded extension operator, and then study u — a that satisfies homogeneous boundary 
conditions. This shows that we can assume that a = 0, then the boundary condition is 
equivalent to w G Hq(M). 

Note first that ||w||hi(m) and ||diu||.L2(M) are equivalent norms on Hq(M), by ( TTBl and the 
Poincare type of inequality for dw, see the first inequality in (fT2|) . The existence part of the 
theorem in the case a = then follows as in [[19, Proposition 5.1.1]. 

To prove the stability estimate, given w G Hq(M), integrate by parts to get 

||u>||^i/C < ||dw||^ 2 = —(Sdw,w) < \\5dw\\ H -i\\w\\ H i. 
That implies (1161) when a = 0, and completes the proof of the lemma. □ 

Remark 1. Lemma [T] in particular justifies the solenoidal-potential decomposition of tensors 
/ with L 2 only regularity, see also [TB]. Then / = f s + dv with v G solving Sdv = Sf. 

3.2. Proof of Theorem [T](b) and Corollary [H To prove Corollary [J we define the set 
Q of simple metrics as follows: near any real analytic simple g, we choose a small enough 
neighborhood in the C 3 (M) topology, so that (j4j) still holds. Clearly, this set is open and 
dense. It remains to prove that this can be done, which is the statement of Theorem (H[b). 

We will prove a bit more. Fix a simple go G C 3 (M) (not necessarily real analytic) and 
assume that (J3]) holds; in particular, the corresponding ray transform I is s-injective. We 
will show that there exists < e 1 so that for any other g with \\g — go\\c 3 (M) < e > © still 
holds with possibly a different constant C > 0, independent on g and e. 

We will apply first [5 S Proposition 4]. There, the weighted ray transform 

W(7)= / w(7(*),700)/(7(s))da 

of functions is studied. The estimate in [5], Proposition 4] compares two such transforms with 
different weights and different metrics (actually, we study more general families of curves 
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in [5]). In our context, we apply Proposition 4] to each fijYi^ before summing up by 
treating 7*7-' as a weight. All we need to show is that the generators of the geodesic flows 
related to go and g are 0(e) close in C 2 . This follows from our assumption \\g — go\\c 3 (M) < £■ 
Then we get 

(18) \\(N 9 -N go )f\\ HHMi) <Ce\\f\\ L , {M) . 

To perturb the l.h.s. of (HI), we need to compare the solenoidal projections f go and f g of / 
related to go and g, respectively. Recall that f s = f — d(6d)]j 5f, where (5d)o is the Dirichlet 
realization of the elliptic operator Sd. By JTBJ, Lemma 1], (5d)~^ : if _1 (M) — > H^(M) 
depends continuously on g G C (M). The same is true for d and S in the corresponding 
spaces. Therefore, 

( 19 ) \\f s 9 -f s J L > m ^ c ° e \\fh*(M)- 

It is enough to prove (@J for / solenoidal, w.r.t. g, so, let us assume that. Then by fTLSj) . and 
our assumption that (@J holds for g , 

Wfgo \\l 2 (m) < C\\N go f\\ HHMl) < Ce\\f\\ LH M) + C\\N g f\\ HHMl) . 

Apply <HM , where f g =fto get 

(1 -C e)||/|| i2(M) < Ce\\f\\ L , (M) + C\\N g f\\ 

Therefore, if e <C 1, we still have (j4j). This completes the proof of Theorem [T](b). Now, 
Corollary [1] follows immediately. 
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